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Abstract 

We  observe  the  quantum  Zeno  effect — where  the  act  of  measurement  slows  the  rate  of  quantum  state 
transitions — in  a  superconducting  qubit  using  linear  circuit  quantum  electrodynamics  readout  and  a 
near-quantum-limited  following  amplifier.  Under  simultaneous  strong  measurement  and  qubit 
drive,  the  qubit  undergoes  a  series  of  quantum  jumps  between  states.  These  jumps  are  visible  in  the 
experimental  measurement  record  and  are  analyzed  using  maximum  likelihood  estimation  to 
determine  qubit  transition  rates.  The  observed  rates  agree  with  both  analytical  predictions  and 
numerical  simulations.  The  analysis  methods  are  suitable  for  processing  general  noisy  random 
telegraph  signals. 


1.  Introduction 

The  backaction  of  measurement  is  a  peculiarly  quantum  mechanical  phenomenon  which  gives  rise  to  striking 
outcomes,  such  as  the  quantum  Zeno  effect  (QZE).  In  the  QZE,  the  act  of  measurement  inhibits  transitions 
between  eigenstates  of  the  measured  observable,  slowing  the  state  evolution  of  a  watched’  quantum  system.  The 
QZE  was  described  in  its  modern  form  in  1977  by  Misra  and  Sudarshan  [1],  although  some  related  questions  were 
tackled  in  prior  papers  [2, 3],  and  the  potential  for  repeated  measurements  to  influence  the  state  evolution  of  a 
quantum  system  was  already  known  to  von  Neumann  in  1932  [4].  The  slowing  of  state  evolution  due  to  the  QZE 
disappears  in  the  classical  limit  H  — >  0,  making  the  QZE  a  useful  test  for  quantum  behavior  in  a  system  [5, 6] . 

The  QZE  was  first  observed  experimentally  in  an  ensemble  of  trapped  ions  [7],  and  has  since  been  seen  in  a 
variety  of  other  systems,  including  the  electronic,  nuclear,  or  motional  states  of  atoms  and  molecules  [8-11], 
optical  photons  [12—14],  microwave  photons  [15],  and  NV  centers  [16].  In  driven  superconducting  qubits,  the 
QZE  has  been  indirectly  inferred  from  the  transition  between  coherent  Rabi  oscillations  and  incoherent 
exponential  population  decay  with  increasing  measurement  strength  [17],  and  by  studying  the  dependence  of 
this  exponential  decay  on  the  time  between  discrete  qubit  projection  pulses  [18].  However,  theoretical  proposals 
also  exist  to  observe  the  QZE  in  the  quantum  trajectory  of  a  continuously  monitored  superconducting  qubit  [19] 
or  in  the  suppression  by  measurement  of  low-frequency  superconducting  flux  qubit  dephasing  [20]. 

The  QZE  and  related  phenomena  can  have  practical  applications  in  quantum  control  and  the  engineering  of 
decoherence.  Carefully  designed  measurements  can  be  used  to  divide  a  larger  Hilbert  space  into  separate  £Zeno 
subspaces’  [21, 22],  where  state  evolution  between  subspaces  is  inhibited  by  the  measurements.  The  resulting 
'quantum  Zeno  dynamics’  have  recently  been  demonstrated  experimentally  [23-25].  This  technique  can  even  be 
used  to  generate  multiparticle  entanglement  directly  [26-29]. 
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In  this  work,  we  report  the  direct  observation  of  the  QZE  in  a  superconducting  qubit  undergoing  continuous 
strong  measurement  with  simultaneous  qubit  driving.  The  measurement  record  is  analyzed  to  extract  quantum 
jumps  indicating  individual  qubit  state  transitions,  and  to  determine  the  rates  at  which  they  occur  in  the 
presence  of  simultaneous  qubit  drive  and  measurement.  The  extracted  transition  rates  show  inhibition  of  qubit 
state  transitions  due  to  the  measurement,  in  agreement  with  both  analytical  Zeno  theory  and  numerical 
simulations.  Additionally,  we  examine  the  excited  state  decay  of  the  qubit  during  measurement,  finding 
qualitative  agreement  with  predictions  for  the  Purcell  decay  of  a  multi-level  qubit  [30, 31].  The  method  for 
extracting  transition  rates  can  also  be  used  to  analyze  general  noisy  random  telegraph  signals. 


2.  Theory  of  the  Zeno  effect 


In  the  absence  of  measurement,  a  resonantly  driven  qubit  will  undergo  sinusoidal  Rabi  oscillations  between 
states  at  frequency  0/27 r,  where  O  depends  on  the  strength  of  the  resonant  drive.  Repeated  projective 
measurements  made  on  this  system  at  time  intervals  r<  I/O  will  tend  to  pin  the  qubit  in  one  state,  with 
occasional  state  transitions  occurring  as  sudden  quantum  jumps.  The  probability  per  unit  time  of  a  quantum 
jump  out  of  the  current  state  is  given  by  (see  e.g.  [32]): 


r  = 


(i) 


where  fm  =  1/r  is  the  frequency  of  the  measurements.  A  similar  result  can  be  derived  in  the  case  of  a  resonantly 
driven  qubit  undergoing  continuous  measurement,  parameterized  by  the  measurement  rate  rm  at  which 
information  is  extracted  from  the  qubit  [33]: 


r  = 


02 

rm' 


(2) 


These  expressions  show  that  the  rate  of  quantum  jumps  between  states  for  a  qubit  being  measured  continuously 
at  rate  rm  is  the  same  as  that  for  a  qubit  undergoing  discrete  measurements  at  time  intervals  r  when  Tmr  =  4. 
This  prediction,  due  to  Schulman  [33],  was  first  demonstrated  experimentally  by  Streed  etal  [10]. 

The  two  major  hallmarks  of  the  QZE  can  be  seen  in  the  expressions  (1)  and  (2).  First,  T  is  independent  of 
time,  meaning  that  the  state  evolution  under  repeated  measurement  is  exponential  (linear  in  time  at  short 
times),  whereas  in  the  absence  of  measurement  the  qubit  would  exhibit  sinusoidal  state  evolution  (quadratic  in 
time  at  short  times).  Qubit  evolution  processes  which  were  already  exponential  in  time,  for  example  decay, 
remain  exponential  and  are  unaffected  by  the  presence  of  repeated  measurements.  Second,  the  transition  rate 
varies  inversely  with  the  measurement  frequency,  tending  to  zero  in  the  limit  of  infinitely  frequent 
measurements.  However,  the  complete  ‘freezing  of  qubit  evolution  does  not  occur  in  real  physical  systems.  This 
is  due  to  the  energy- time  uncertainty  relationship,  which  causes  the  qubit  to  couple  to  an  arbitrarily  large  energy 
spectrum  (and  thus  arbitrarily  many  decay  channels)  as  the  time  between  measurements  goes  to  zero  [34, 35]. 

Our  experiment  uses  a  superconducting  qubit  in  the  circuit  QED  architecture  [36].  We  bias  the  qubit  in  the 
dispersive  regime,  where  the  qubit-cavity  detuning  A  is  much  larger  than  the  qubit-cavity  coupling  g.  In  this 
limit,  the  cavity  resonance  frequency  depends  on  the  qubit  state,  enabling  a  quantum  non- demolition 
measurement  of  the  qubit  state  by  driving  the  cavity  near  resonance  and  observing  the  response  of  its  steady- 
state  field.  The  measurement  strength  is  directly  related  to  the  distinguishability  of  the  different  cavity  field 
amplitudes  conditioned  on  the  qubit  states.  For  a  linear  cavity,  the  distinguishability  is  proportional  to  the 
amplitude  of  the  measurement  drive  and  could  in  principle  be  increased  ad  infinitum.  In  real  systems,  however, 
cavity  nonlinearities  (both  intrinsic  and  induced  by  coupling  to  the  anharmonic  qubit  circuit)  limit  the 
obtainable  field  separations  for  readout,  even  at  moderate  drive  amplitudes  [37]. 

Following  [19],  and  as  shown  in  appendix  A,  we  derive  the  transition  rate  of  a  two-level  qubit  coupled  to  a 
cavity  in  the  Zeno  regime  starting  with  the  Jaynes-Cummings  Hamiltonian  in  the  dispersive  regime  with  an 
additional  qubit  drive: 


H  =  cJrO^a  H - (a;q  +  x)°z  +  X^a(7z  H - (eroaPelujTot  +  6*0ae  luJrot )  H - (£2cr_e  luJdt  +  Q*cr+eluJit). 


(3) 


Here,  a  (aQ  is  the  annihilation  (creation)  operator  for  photons  in  the  cavity,  the  qubit  is  described  by  the  Pauli 
matrices  cr,  and  we  have  made  a  rotating  wave  approximation  in  the  two  driving  terms.  The  parameters  are  the 
cavity  resonance  frequency  o;r,  the  qubit  frequency  cjq,  the  qubit-cavity  coupling  g,  and  the  dispersive  shift 
X  =  g2  /  A,  where  A  =  ujy  —  uq.  We  apply  a  measurement  drive  to  the  cavity  with  frequency  cjro  and  strength 
6ro  and  a  Rabi  drive  to  the  qubit  with  frequency  and  amplitude  0.  We  then  perform  a  polaron  transformation 

on  the  system,  which  effectively  describes  the  cavity  in  a  semi- classical  picture  with  coherent  state  amplitudes 
that  are  conditioned  on  the  qubit  state.  This  transformation  enables  us  to  decouple  the  dynamical  equations  of 
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Figure  1.  Experimental  setup.  A  schematic  of  the  experimental  setup  is  shown  in  (a).  The  readout  and  qubit  drive  tones  are  sent  to  the 
qubit/ cavity  circuit  via  a  weakly  coupled  input  port.  The  readout  signal  from  the  cavity  exits  through  the  strongly  coupled  port  and  is 
amplified  by  a  superconducting  parametric  amplifier,  as  well  as  further  cryogenic  and  room  temperature  amplifiers  (not  shown), 
before  being  detected  with  homodyne  mixing  and  digitized.  The  Zeno  experiment  pulse  sequence  is  shown  in  (b).  The  readout  is 
turned  on  3  (is  before  the  qubit  drive  to  allow  the  cavity  to  come  to  steady  state.  Only  the  data  from  the  10  (is  with  both  qubit  drive  and 
readout  are  processed  for  Zeno  rate  analysis. 


the  qubit  and  cavity  and  obtain  a  qubit- only  reduced  master  equation,  where  the  qubit  dynamics  depend  mainly 
on  qubit  drive  and  measurement  strength.  Because  of  the  presence  of  the  Rabi  drive,  this  treatment  is  only  valid 
in  the  case  of  small  measurement  strength,  when  the  cavity  field  amplitudes  corresponding  to  the  two  qubit 
states  are  nearly  indistinguishable.  Solving  this  master  equation  yields  a  qubit  transition  rate  from  the  ground 
state  to  excited  state  in  the  presence  of  continuous  circuit  QED  measurement  [19]: 


n2 


G,  drive 


2(72  +  Id) 


(4) 


where  1/  is  the  measurement- induced  dephasing  rate  and  y2  is  the  intrinsic  qubit  dephasing  rate  in  the  absence 
of  measurement.  The  measurement  rate  is  defined  as  rm  =  21^,  as  in  [19].  In  the  strong  measurement  limit 
where  1/  72,  this  expression  is  the  same  as  equation  (2).  Although  equation  (4)  is  derived  in  the  limit  where 

the  cavity  linewidth  k  rm,  the  similarity  to  equation  (2) — which  has  no  such  constraint  in  its  derivation — 

suggests  that  the  result  may  still  be  valid  when  k  <  rm,  which  is  the  regime  of  our  experimental  data. 


3.  Experimental  setup  and  calibrations 

The  experimental  setup  is  shown  schematically  in  figure  1  (a).  Our  qubit  is  a  planar  two-junction  transmon  qubit 
[38],  tuned  with  a  dc  magnetic  flux  to  a  Lamb-shifted  qubit  frequency  of  £uq/27r  =  5.3556  GHz  and  an 
anharmonicity  of  a/2n  =  —  258  MHz.  We  denote  the  two  lowest  transmon  energy  levels  as  |g)  and  |  e), 
respectively.  The  qubit  is  capacitively  coupled  (g/ 2tt  =  105.3  MHz)  to  a  lumped- element  readout  cavity  with 
bare  frequency  of  cjy/2tt  =  6.2724  GHz  and  linewidth  ^/27r  =  7  MHz.  Qubit  measurement  was  performed  by 
applying  a  readout  drive  to  the  cavity  at  ujyo/ 2i r  =  6.282  GHz.  The  readout  signal  from  the  cavity  was  amplified 
with  a  near- quantum-limited  superconducting  parametric  amplifier  [39]  followed  by  additional  cryogenic  and 
room  temperature  amplification  stages,  detected  with  homodyne  mixing,  and  digitized  at  10  ns  intervals.  Both 
the  qubit/ cavity  system  and  the  parametric  amplifier  were  anchored  to  the  mixing  chamber  of  a  dilution 
refrigerator  at  50  mK. 

The  average  cavity  photon  occupation  fi  and  the  measurement- induced  dephasing  rate  I/,  which 
characterize  the  measurement  strength,  were  determined  by  qubit  spectroscopy  with  a  simultaneous  readout 
tone.  The  value  of  fi  was  determined  from  the  ac  Stark  shift  of  the  qubit  frequency,  accounting  for  qubit- induced 
cavity  nonlinearities  [30].  Because  of  the  choice  of  readout  frequency,  fi  is  different  for  the  two  qubit  states  |g) 
and  |  e) ;  the  ac  Stark  shift  measurements  give  the  average  photon  occupation  with  the  qubit  in  the  ground  state, 
ng,  since  the  qubit  spectroscopy  tone  was  much  weaker  than  the  saturation  amplitude  of  the  qubit  transition.  We 
use  numerical  simulations  to  infer  the  value  of  fie  as  a  function  of  the  measured  fig.  The  value  of  1/  is  given  by  the 
half- width  at  half-maximum  (HWHM)  of  the  ac- Stark- shifted  qubit  line  [40],  which  was  determined  by  fitting 
to  the  measured  line  shape  for  several  qubit  drive  powers  and  extrapolating  to  find  the  HWHM  at  zero  qubit 
drive  power.  This  method  allows  us  to  calibrate  out  the  effects  of  power  broadening  due  to  the  qubit  drive  tone 
on  the  measured  qubit  linewidth.  With  the  readout  tone  off  during  qubit  evolution,  we  measured  a  qubit 
relaxation  time  of  2J  =  575  ns  and  a  pure  dephasing  time  of  =  7.9  /is,  corresponding  to  an  intrinsic 

dephasing  rate  72  =  1/2  2J  +  1/2/ =  1  /is-1. 

The  qubit  drive  frequency  was  chosen  to  be  the  ac- Stark- shifted  qubit  frequency  ujq  (fig).  Since  the  qubit 

drive  tone  reaches  the  qubit  via  the  cavity,  its  amplitude  at  the  qubit  (and  thus  the  no-measurement  Rabi 
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Figure  2.  Data  processing  and  quantum  jump  extraction.  Part  (a)  is  a  histogram  of  detected  voltages  from  the  full  data  set  (104  readouts 
with  103  points  each)  with  Tm  =  134  fis~l  (ng  =  3.4)andD/27r  =  3.6  MHz,  both  for  the  raw  data  (red)  and  after  filtering  (black). 
The  hysteretic  threshold  voltages  for  state  determination  are  shown  as  horizontal  dashed  lines.  In  (b),  a  typical  single  measurement 
trace  from  the  same  data  set  is  shown  both  as  raw  data  (red)  and  after  filtering  (black),  with  the  state  determination  threshold  voltages 
shown  as  horizontal  dashed  lines.  The  background  color  indicates  the  extracted  qubit  state:  white  for  |g)  and  blue  for  \e).  Panels  (c) 
and  (d)  show  equivalent  data  for  Tm  =  134  ^s-1  and  D/27 r  =  0.8  MHz.  Because  the  transition  rates  are  lower  than  in  (a)  and  (b),  the 
extraction  algorithm  performs  heavier  filtering,  resulting  in  larger  filtered  SNR  and  smaller  hysteresis  between  the  threshold  voltages. 


frequency  £2)  depends  on  both  its  frequency  and  its  amplitude  at  the  cavity  input.  We  calibrated  £2  by  measuring 
the  frequency  of  Rabi  oscillations  (performed  with  measurement  off  during  the  qubit  drive)  as  a  function  of 
qubit  drive  amplitude  for  a  given  qubit  frequency.  We  performed  this  calibration  at  several  different  qubit 
frequencies  spanning  the  range  of  the  lj &  values  used  in  the  Zeno  experiments  by  tuning  the  qubit  frequency  with 
an  external  flux  bias.  We  interpolated  between  these  data  points  to  determine  the  value  of  £2  for  a  given  uj&  and 
qubit  drive  amplitude. 

Each  iteration  of  the  experiment  consisted  of  a  single  measurement  pulse  lasting  17.5  /is,  with  simultaneous 
qubit  drive  applied  for  10  /is  beginning  3  fis  after  the  measurement  was  turned  on,  as  shown  in  figure  1(b).  These 
times  are  chosen  to  be  much  longer  than  the  relevant  timescales  1/ft,  1  /Tm,  l/£2,  and  T),  so  that  the 
experimental  record  captures  the  steady- state  dynamics  of  the  system.  We  recorded  104  such  iterations,  spaced 
100  /is  apart,  for  each  combination  of  measurement  strength  and  qubit  drive  strength.  The  measurement  record 
was  digitized  at  108  samples  s-1. 

4.  Data  processing  and  analysis 

The  measurement  records  described  above  were  analyzed  to  extract  the  transition  rates  from  | g)  — >  |  e),  denoted 
Ip  and  from  |  e)  — >  |g),  denoted  I}.  With  multiple  gigabytes  of  raw  quantum  jump  data  to  analyze,  the  data 
processing  algorithm  must  be  able  to  operate  with  minimal  user  input  and  provide  reliable  output  over  a  broad 
range  of  signal-to-noise  ratio  (SNR)  and  qubit  transition  rates,  including  low- SNR,  high-rate  scenarios. 

We  first  perform  filtering  of  the  raw  data  to  reduce  the  noise  bandwidth  and  increase  the  SNR.  We  use  a  zero- 
delay  Gaussian  finite-impulse-response  filter  to  ensure  that  state  transition  edges  will  remain  smooth  and  will 
not  be  shifted  in  time  by  the  filtering  process.  The  order  of  the  filter  must  be  chosen  to  reduce  the  noise  without 
degrading  the  desired  signal,  and  so  will  depend  both  on  the  initial  SNR  and  the  transition  rate  between  states. 
For  each  experimental  bias  point,  we  histogram  all  107  data  points  taken  during  simultaneous  qubit  driving  and 
readout  (104  traces  of  103  points  each).  Example  histograms  of  raw  data  are  shown  as  the  red  curves  in 
figures  2(a)  and  (c).  We  then  perform  filtering  and  histogram  the  filtered  data  points,  shown  as  black  curves  in 
figures  2(a)  and  (c).  The  histogram  of  the  filtered  data  forms  a  bimodal  distribution  ,  with  a  minimum  between 
the  two  peaks  corresponding  to  the  two  qubit  states.  The  height  of  this  minimum  is  dictated  both  by  the  amount 
of  noise  on  the  signal  and  the  number  of  transitions  in  the  measurement  record.  Filtering  will  reduce  the  noise 
amplitude,  and  thus  the  width  and  overlap  of  the  bimodal  peaks,  leading  to  a  reduction  in  the  height  of  the 
minimum.  However,  filtering  will  also  slow  the  rise  and  fall  times  of  state  transition  edges,  which  will  increase 
the  height  of  the  minimum  between  histogram  peaks.  Based  on  these  considerations,  we  designate  the  ‘optimal’ 
filter  order  as  the  one  which  gives  the  lowest  minimum  between  the  bimodal  peaks  of  the  resulting  voltage 
histogram,  and  determine  it  by  applying  filters  of  increasing  order  to  the  raw  data  until  the  height  of  the 
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If  the  histogram  is  not  bimodal,  as  for  the  raw  data  in  figure  2(c),  we  perform  additional  filtering  until  it  becomes  bimodal. 
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minimum  between  the  peaks  is  as  small  as  possible.  This  provides  a  simple,  robust  method  for  choosing  filter 
order.  The  ‘optimal’  filter  order  was  determined  separately  for  each  combination  of  qubit  drive  strength  and 
measurement  strength. 

After  the  data  have  been  filtered,  we  determine  the  qubit  state  at  each  time  point  using  a  thresholding 
algorithm.  For  increased  robustness  to  noise,  we  use  two  hysteretic  thresholds  with  ‘Schmitt  trigger’  behavior 
[4 1  ];  state  transitions  are  registered  when  the  voltage  crosses  the  higher  threshold  going  upwards  (if  the  state  was 
low),  or  when  it  crosses  the  lower  threshold  going  downwards  (if  the  state  was  high).  The  voltage  thresholds  are 
shown  as  dotted  horizontal  lines  in  figure  2.  Figures  2(b)  and  (d)  show  representative  individual  data  traces  from 
the  data  sets  histogrammed  in  figures  2(a)  and  (c),  respectively,  corresponding  to  two  different  values  of  D  with 
identical  Tm.  The  red  curves  show  the  raw  data,  while  the  black  curves  are  after  filtering.  The  blue  and  white 
background  colors  show  the  extracted  qubit  state  (|  e)  and  \g),  respectively)  as  determined  by  the  thresholding 
algorithm.  The  data  in  (c)  and  (d)  have  fewer  state  transitions,  and  thus  heavier  filtering  can  be  employed  to 
increase  the  SNR. 

The  location  of  the  thresholds  is  determined  using  histograms  of  the  filtered  voltages,  and  depends  on  the 
voltage  V"m  of  the  minimum  between  the  histogram  peaks  and  the  voltages  Vh  ( VJ)  and  the  HWHM  wh  (wi)  of  the 
high  (low)  histogram  peaks.  The  value  of  w h  (wf)  is  derived  from  the  half-height  of  the  histogram  above  (below) 
the  corresponding  peak  yh  ( Vj),  a  method  that  works  even  when  the  histogram  peaks  are  not  well- separated.  We 

define  the  SNR  of  the  filtered  data  to  be  y/2  In  2  ( —  Vi)  / (wh  +  w\).  The  thresholds  are  chosen  to  be 

2  2 

Vm  +  2ln2™y — and  Vm  —  2]n2^ — — ,  as  suggested  in  [4 1  ] ,  although  unlike  that  work  we  do  not  perform 
iterative  refinement  of  the  thresholds.  The  distance  between  the  thresholds  increases  as  the  SNR  decreases  to 
reduce  the  likelihood  of  spurious  state  transitions  being  registered  due  to  noise.  For  each  measurement  trace,  we 
determine  the  initial  state  by  comparing  the  first  data  point  with  Vm.  Given  a  low  (high)  state,  we  look  for  the  first 
upward  crossing  of  the  higher  threshold  (downward  crossing  of  the  lower  threshold)  and  note  the  dwell  time 
before  the  state  change.  We  continue  in  this  manner,  alternating  thresholds,  until  the  end  of  the  trace  is  reached, 
noting  down  all  dwell  times  in  both  the  high  and  low  states.  For  each  measurement  trace,  exactly  one  dwell  time 
is  cut  short  by  the  end  of  the  trace  and  is  marked  as  ‘right-censored’  [42] . 

We  use  maximum  likelihood  estimation  to  determine  the  transition  rates  between  states  given  the  set  of 
observed  dwell  times  in  each  state.  However,  the  filtering  can  cause  some  fast  voltage  excursions  in  the  raw  data, 
like  the  one  in  figure  2(d)  near  8.25  /is,  not  to  register  as  state  transitions.  In  general,  the  finite  bandwidth  of  the 
measurement  chain,  combined  with  filtering  described  above,  skews  the  distribution  of  observed  dwell  times 
toward  longer  times.  We  compensate  for  this  in  our  analysis  by  assuming  a  probability  distribution  for  the  dwell 
times  which  takes  these  effects  into  account  [43].  For  right- censored  dwell  times,  the  observed  time  represents  a 
lower  bound  on  the  dwell  time,  rather  than  an  exact  value.  However,  we  can  include  this  partial  information  in 
the  likelihood  function  as  well,  which  is  particularly  useful  for  data  sets  with  relatively  few  state  transitions  and 
thus  a  higher  fraction  of  censored  dwell  times.  We  emphasize  that  this  maximum  likelihood  method  can  be 
applied  to  dwell  times  extracted  using  any  technique;  for  example,  it  can  be  used  with  dwell  times  determined 
from  wavelet  analysis,  which  is  more  robust  than  thresholding  for  signals  with  substantial  low-frequency  drifts 
[44].  Details  on  the  maximum  likelihood  estimation  and  functional  forms  of  the  likelihood  function  are 
provided  in  appendix  D. 

The  extraction  algorithm  was  tested  on  simulated  noisy  random  telegraph  data  with  a  variety  of 
experimentally  relevant  transition  rates  and  SNRs.  We  note  that  these  simulated  data  are  not  the  same  as  the 
numerical  simulations  of  the  qubit  dynamics  described  in  the  next  section.  The  3  dB  noise  bandwidth  (14  MHz) 
and  sample  rate  (108  sample  s-1)  were  the  same  as  in  the  raw  experimental  data.  The  error  in  the  extracted  mean 
transition  rate  (F)  +  T|)  /2,  relative  to  the  true  mean  transition  rate,  is  plotted  in  figure  3  for  four  values  of  the 
initial  SNR  before  filtering.  For  simulated  data  at  the  highest  rates  shown,  the  mean  dwell  time  in  a  given  state 
was  130  ns,  while  the  shortest  oscillation  timescale  of  the  simulated  noise  was  around  70  ns,  making  the  task  of 
distinguishing  signal  from  noise  challenging.  Despite  this,  the  mean  systematic  bias  in  the  extracted  rates  is 
below  12%,  even  at  an  SNR  of  only  1.5.  This  performance  is  also  notable  because  the  algorithm  operates 
autonomously,  without  the  need  for  externally  provided  guesses  or  input,  across  the  entire  range  of  transition 
rates  and  SNR  shown. 

The  transition  rates  extracted  from  the  simulated  noisy  random  telegraph  data  were  used  to  calibrate 
systematic  bias  and  systematic  uncertainty  in  the  rates  extracted  from  the  experimental  data,  by  matching 
experimental  traces  to  simulated  data  sets  with  similar  extracted  transition  rates  and  filtered  SNR.  The  median 
systematic  bias  and  median  systematic  uncertainty  for  the  experimental  data  were  —3.5%  and  2%  of  the 
extracted  rates,  respectively.  The  magnitude  of  the  systematic  bias  (uncertainty)  was  below  20%  (10%)  of  the 
extracted  rate  for  all  experimental  data  points,  and  below  9%  (5%)  of  the  extracted  rate  for  90%  of  the 
experimental  data  points. 
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5.  Numerical  simulations 

We  also  performed  numerical  simulations  of  the  qubit/ cavity  system  over  a  range  of  bias  points,  using  the 
generalized  Jaynes-Cummings  Model  for  a  multi-level  qubit  coupled  to  a  cavity.  We  simulated  the  dynamics 
numerically,  including  the  measurement,  with  a  stochastic  master  equation  [19, 45, 46],  using  experimentally 
measured  parameters  for  the  qubit  and  cavity.  The  master  equation  included  qubit  relaxation  and  dephasing  as 
well  as  cavity  photon  loss.  The  measurement  and  qubit  drive  were  simulated  as  two  independent  coherent 
microwave  drives  acting  on  the  cavity.  The  simulations  yield  measurement  records  for  direct  comparison  with 
experiment,  as  well  as  records  of  qubit  populations,  complex  cavity  field  amplitudes,  and  cavity  photon 
occupation  numbers.  Details  are  given  in  appendix  B.  Qubit  transition  rates  were  then  determined  from  the 
simulated  qubit  state  population  record  using  the  same  rate  extraction  algorithm  used  to  process  the 
experimental  data. 

The  experimental  data,  and  thus  the  rate  extraction  algorithm,  do  not  distinguish  between  the  qubit  state  |  e) 
and  higher  excited  states.  This  occurs  because  we  operate  the  parametric  amplifier  in  phase- sensitive  mode  to 
achieve  the  lowest  noise  performance,  and  thus  only  one  quadrature  of  the  readout  signal  from  the  cavity  is 
amplified  [39] .  The  phase  of  the  amplified  quadrature  was  chosen  to  maximize  the  ability  to  discriminate 
between  states  |g)  and  |  e) ;  however,  because  of  the  parameters  of  the  qubit/ cavity  system,  the  projection  of  the 
complex  cavity  amplitudes  corresponding  to  higher  excited  states  onto  this  choice  of  amplified  quadrature  was 
essentially  the  same  as  that  for  |  e),  rendering  them  indistinguishable  from  |  e)  in  the  measurement  record  [47]. 
The  populations  of  higher  excited  states  were  measured  to  be  ^  1%  —  5%  using  a  different  data  set  which 
distinguished  |  e)  from  higher  qubit  states  by  choosing  a  different  amplified  quadrature,  at  the  expense  of 
substantially  decreased  SNR.  These  values  are  corroborated  by  the  numerical  simulations.  However,  the 
presence  of  qubit  population  outside  the  { \g) ,  |  e) }  manifold  is  not  expected  to  have  an  effect  on  the  estimates 
ofrT. 

6.  Results 

We  first  examine  the  qubit  excitation  rate  during  simultaneous  measurement  and  qubit  driving  to  look  for  the 
QZE.  Since  fig  fie  for  our  experimental  parameters,  the  qubit  drive  tone  at  the  ac- Stark- shifted  qubit 

frequency  ujq  (ng)  was  resonant  when  the  qubit  was  in  the  ground  state,  but  not  when  it  was  in  the  excited  state. 
For  this  reason,  our  analysis  of  driven  transition  rates  is  restricted  to  the  |g)  — >  |  e)  transition  rate,  denoted  If, 
which  can  be  expressed  as  the  sum  of  three  rates: 

If  —  If, drive  +  I^DD  +  I^th-  (5) 

Here  Thrive  is  the  transition  rate  due  to  the  qubit  drive,  Tf  >DD  is  the  contribution  from  dressed  dephasing 
[37, 48],  and  F/  th  represents  thermal  excitation  of  the  qubit.  To  isolate  Thrive  f°r  comparison  to  the  predictions 
in  equation  (4),  we  determine  1}  as  a  function  of  measurement  strength  in  the  absence  of  qubit  drive  (D  =  0), 
which  is  equal  to  Tf  jDD  +  If  jth.  We  can  then  subtract  this  contribution  from  the  total  rate  If  in  the  presence  of 
qubit  drive  to  find  If  5 drive  -  This  calibration  was  performed  for  both  experimental  data  and  numerically  simulated 
data.  The  value  of  Tf  jDD  +  If  jth  was  measured  to  be  between  0.018  and  0.022  /is-1  for  all  bias  points,  which  is 
considerably  smaller  than  If  5  drive  f°r  all  but  the  lowest  values  of  Q. 

Figures  4  and  5  show  the  extracted  Thrive as  a  function  of  D  for  ten  different  measurement  strengths  ranging 
from  Tm  =  134  /is-1  to  Tm  =  393  /is-1,  corresponding  to  values  of  ng  between  3.4  and  37.  The  rates  from  the 
experimental  data  and  the  numerical  simulations  are  plotted  as  red  squares  and  black  circles,  respectively,  while 
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Figure  4.  Driven  transition  rates  during  measurement.  We  plot  values  of  FV  drive  versus 

qubit  drive  strength  O  for  ten  measurement 

strengths  Tm  ranging  between  134  /i s-1  (ng  =  3.4)  and  393  fis~l  ( ng  =  37).  Experimental  data  (red  squares)  are  shown  along  with 
numerically  simulated  data  (black  circles)  and  theoretical  values  calculated  from  equation  (4)  with  no  adjustable  parameters  (blue 
lines).  Error  bars  represent  95%  confidence  intervals.  Horizontal  error  bars  are  smaller  than  the  symbol  widths  and  are  not  shown  for 

clarity. 

(a)  (b) 


(c) 


(h) 


(d) 


(o) 


(k) 


Figure  5.  Driven  transition  rates  during  measurement  on  logarithmic  axes.  The  data  are  the  same  as  in  figure  4,  but  the  logarithmic 
axes  highlight  the  agreement  between  experiment,  analytical  theory,  and  numerics  for  the  smallest  Q,  as  well  as  the  Q?  scaling  expected 
from  equation  (4).  Experimental  data  (red  squares)  are  shown  along  with  numerically  simulated  data  (black  circles)  and  theoretical 
values  calculated  from  equation  (4)  with  no  adjustable  parameters  (blue  lines).  Error  bars  represent  95%  confidence  intervals. 
Horizontal  error  bars  are  smaller  than  the  symbol  widths  and  are  not  shown  for  clarity. 


the  rates  predicted  by  the  analytical  theory  in  equation  (4),  with  no  adjustable  parameters,  are  plotted  as  solid 
blue  lines. 

The  experimental  data  and  numerics  are  in  good  agreement,  and  both  coincide  with  the  prediction  of  the 
analytical  theory  at  most  points,  showing  the  presence  of  the  QZE  across  a  broad  range  of  measurement 
strengths  and  qubit  drive  amplitudes.  Some  deviations  can  be  seen  above  Vt/2i:  «  6  MHz,  where  both  the 
experiment  and  numerics  give  consistently  lower  values  of  If  ? drive  than  predicted  by  the  analytical  theory. 
Because  the  numerical  simulations  exhibit  the  same  behavior  as  the  experiment,  we  believe  this  is  a  real  effect 
and  not  an  experimental  artifact.  We  postulate  that  this  additional  slowing  of  the  transition  rate  maybe  due  to 
higher  order  terms  neglected  in  the  derivation  of  the  analytical  theory  under  the  approximation  Tm  <C  k  ,  since 
this  approximation  is  no  longer  valid  here.  The  effect  is  somewhat  more  pronounced  for  larger  Tm.  This  could 
also  be  associated  with  the  breakdown  of  the  dispersive  approximation  for  ng  >  nCY[t  =  A2/4g2,  which 
corresponds  in  our  system  to  Tm  Pd  340  fi s-1. 

We  additionally  examine  the  qubit  Tx —  the  time  constant  for  transitions  from  |  e)  to  | g) — with 
measurement  on  but  no  qubit  drive.  In  the  simple  two-level  picture  used  to  derive  equation  (4),  transition  rates 
from  exponential  processes  such  as  decay  are  unchanged  by  the  QZE.  However,  in  our  driven  multilevel 
system  other  effects  are  expected  to  give  some  dependence  of  Tx  on  measurement  strength.  In  circuit  QED,  the 
T i  decay  rate  can  in  general  be  parameterized  as  the  sum  of  the  Purcell  decay  rate  TP  and  a  decay  rate  TNR  from 
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Figure  6.  Tj  versus  ne  with  no  qubit  drive.  The  red  squares  are  measured  values  of  the  qubit  decay  rate  as  a  function  of  ne.  The  black 
triangles  and  circles  are  values  derived  from  numerical  simulations;  the  former  distinguishes  between  state  |  e)  and  higher  states  in  the 
rate  extraction  process,  while  the  latter  does  not  (as  is  the  case  for  our  experimental  data).  These  values  differ  somewhat  for  ne  >  8, 
indicating  the  presence  of  population  in  higher  excited  states.  The  solid  red  line  is  the  predicted  Purcell  decay  rate  TP  from  equation  (6), 
and  the  dashed  red  line  is  Tp  +  E^r  as  described  in  the  text.  The  blue  circles  show  scaling  of  the  qubit  decay  rate  that  would  be 
expected  from  equation  (4)  if  the  decay  were  subject  to  the  QZE.  Error  bars  represent  95%  confidence  intervals. 


nonradiative  loss  channels  such  as  dielectric  loss  [49, 50].  The  Purcell  effect  can  be  thought  of  as  the  decay  of  the 
photonic  part  of  the  qubit  eigenstate  in  the  coupled  qubit-cavity  system,  as  seen  in  the  expression  for  the  Purcell 
rate  [30, 31] 


rP  =  k  |  (g,  n\a  \e,  n)  |2,  (6) 


where  \g,  n)  and  |  e,  n)  are  eigenstates  of  the  generalized  Jaynes-Cummings  Hamiltonian  for  a  multilevel  qubit 
corresponding  most  closely  to  n  photons  in  the  bare  cavity  and  the  bare  qubit  in  |g)  or  |  e),  respectively.  The 
matrix  elements  in  equation  (6)  can  be  calculated  numerically,  and  become  smaller  for  increasing  n.  As  a  result, 
the  Purcell  decay  of  the  qubit  is  suppressed  by  the  presence  of  photons  in  the  readout  cavity  [30,31].  Note  that  in 
our  case  we  replace  n  with  ne,  because  Purcell  decay  from  |  e)  occurs  with  ne  photons  in  the  cavity  initially. 
Details  of  the  calculation  are  given  in  appendix  C. 

Figure  6  shows  the  qubit  decay  rate  as  a  function  of  ne.  The  expected  decay  rate  TP  due  to  the  Purcell  effect,  as 
calculated  from  equation  (6),  is  plotted  as  a  red  solid  line.  The  measured  7\  decay  of  the  qubit  (red  squares), 
which  shows  some  suppression  with  increasing  ne,  is  larger  than  Tp,  indicating  that  non-radiative  loss  channels 
are  contributing  to  the  qubit  decay  (i.e.  TNR  >  0).  We  derive  a  value  for  TNR,  which  we  assume  to  be 
independent  of  ne,  by  subtracting  the  zero-photon  Purcell  decay  rate  Kg2 /  A2  from  the  experimentally 
determined  value  of  =  1  /IJ  for  ne  =  0.  The  theoretically  predicted  total  decay  rate  TP  +  TNR  is  plotted  as  a 
red  dotted  line.  This  prediction  is  reasonably  close  to  the  experimental  data  for  small  ne,  but  deviates 
substantially  with  increasing  ne.  To  determine  whether  this  is  an  experimental  artifact  (due  to  readout-power- 
dependent  dielectric  loss  in  the  transmon  capacitor  [51],  for  example,  which  would  give  an  nt  dependence  to 
TNR),  we  also  performed  numerical  simulations  of  the  qubit  decay  with  no  Rabi  drive.  The  resulting  rates  (black 
circles  and  black  triangles)  are  also  consistently  lower  than  the  theoretical  prediction  for  rP  +  rNR,  and  are  in 
reasonable  agreement  with  the  experimental  data.  Because  the  numerical  calculations  assume  that  TNR  is 
independent  of  ne,  this  provides  confirmation  that  the  suppression  in  Tx  decay  is  not  an  experimental  artifact. 
The  numerics  also  allow  us  to  determine  the  effect  of  qubit  population  in  states  higher  than  |  e)  on  the  observed 
decay  rate;  the  black  circles  show  the  observed  7\  decay  rate  when  |  e)  is  not  distinguished  from  higher  qubit 
states,  as  is  the  case  with  our  experimental  data,  while  the  black  triangles  show  the  true  Tx  decay  rate  from 
I e)  ~ >  I g)  when  |  e)  is  distinguished  from  higher  qubit  states.  The  difference  between  these  traces  becomes 
significant  for  «e>  8  ,  and  shows  that  the  inability  to  distinguish  higher  qubit  states  from  |  e)  in  this  system 
causes  the  qubit  decay  rate  to  appear  lower  than  its  true  value,  but  the  effect  does  not  appear  to  be  large  enough  to 
fully  account  for  the  discrepancy  between  the  experimental  data  and  the  analytical  theory.  We  believe  that  other 
effects,  such  as  dressed  dephasing  [37, 48]  and  qubit- induced  nonlinearities,  are  responsible  for  this  deviation 
from  the  simple  theory  in  equation  (6).  The  blue  circles  show  how  the  qubit  decay  rate  would  depend  on  ne  if 

it  were  subject  to  the  QZE  and  scaled  as  1  /Tm,  as  in  equation  (4),  relative  to  the  measured  value  for  the  lowest  ne. 
The  large  disparity  between  these  points  and  the  experimental  and  numerical  data  indicate  that  the  variation  of 
Ij  with  ne  is  not  due  to  the  QZE;  this  agrees  with  the  expected  result  that  the  QZE  does  not  affect  exponential 
processes  such  as  Tx  decay. 
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7.  Conclusions 

We  have  observed  the  QZE  in  a  continuously  measured  superconducting  qubit  and  demonstrated  quantitative 
agreement  with  both  analytical  theory  and  numerical  simulations.  Interestingly,  the  agreement  holds  even  for 
measurement  strengths  much  larger  than  are  allowed  by  the  assumptions  of  the  analytical  derivation.  We  have 
shown  that  7\  decay  is  not  subject  to  the  QZE,  as  predicted.  We  do  observe  an  gradual  increase  in  the  qubit  7\ 
with  increasing  measurement  strength,  in  qualitative  agreement  with  the  simple  picture  from  the  Purcell  effect 
in  circuit  QED.  We  have  also  demonstrated  a  robust  algorithm  for  determining  transition  rates  from  our 
measurement  record,  which  can  be  readily  applied  to  the  analysis  of  general  noisy  random  telegraph  signals. 
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Appendix  A.  Details  of  analytical  Zeno  derivation 


Here  we  derive  the  expression  for  qubit  transition  rates  given  in  equation  (4).  We  will  focus  only  on  the  steps 
necessary  to  derive  the  qubit  transition  rates;  more  detail  can  be  found  in  [19].  We  note  that  this  derivation  is  for 
the  case  of  a  two-level  qubit,  not  a  multi-level  qubit  such  as  a  transmon.  We  describe  the  qubit-cavity  system  in 
the  Jaynes-Cummings  model  (setting  h  =  1): 

Hjc  =  uT($a  +  uqaz  +  g(<$o-  +  aa+)  +  ~(6™  e-la;rotat  +  6*elcJrota),  (Al) 


where  ljy  is  the  cavity  frequency,  ujq  the  qubit  energy,  g  is  the  qubit-cavity  coupling,  and  ero  is  the  amplitude  of 
the  cavity  drive.  Here  a  is  the  annihilation  operator  for  cavity  photons,  and  the  Pauli  matrices  a  describe  the 
qubit.  For  large  qubit-cavity  detuning,  such  that  A  =  cjy  —  uq  g,  we  move  into  the  dispersive  frame  [36] 


Hdisp  =  (u;r  -  cjro)afa  +  ^Cjq az  +  xaWz  +  |(6roa!  +  4/0, 


(A2) 


where  x  =  g2/  A  is  the  dispersive  shift  and  Cbq  =  cjq  +  x  is  the  Lamb-shifted  qubit  frequency.  We  also  move 
the  cavity  into  a  rotating  frame  at  frequency  cjro  and  assume  that  cjro  is  far  detuned  from  the  qubit  transition  ujq. 
For  a  cavity  drive  used  for  a  dispersive  measurement,  this  condition  is  fulfilled.  Since  we  are  concerned  with  the 
qubit’s  dynamics  under  simultaneous  qubit  drive  and  measurement,  we  add  a  qubit  driving  term  at  frequency 
to  equation  (A2),  writing 

Hqi  =  +  f 2*o-+eiw“f),  (A3) 


where  Q  is  the  qubit  Rabi  frequency.  Hereafter  we  examine  the  case  of  resonant  qubit  drive,  where  =  Coq.  In  a 

frame  rotating  at  the  qubit  drive  frequency,  the  uncoupled  qubit  dynamics  are  described  by  Hq  =  leading 

to  coherent  qubit  flopping  at  the  Rabi  frequency.  The  complete  dynamics  of  the  system  will  be  described  by  a 
master  equation  of  the  form 

P  =  —  i [TYdisp  +  #qd>  p\  +  KV[a\p  +  7i V[a-]p  +  j ipV[crz\p, 


with  the  cavity  loss- rate  k,  the  qubit  decay  rate  and  the  qubit  pure  dephasing  rate  7  .  The  total  qubit 

dephasing  rate  is  then  y2  =  ^  +  7^.  The  dissipative  superoperators  are  defined  by 

V[d]p  =dpd^  -  ^(d^dp  +  pd^d). 

Next  we  apply  a  polaron  transformation  to  the  system,  defined  by  [  1 9] 

P  =  YJD(ui)  |i)(i|,  (A4) 

i 

where  D  (ce)  =  exp  { acfi  —  a* a]  is  the  usual  field  displacement  operator  and  the  sum  is  over  the  qubit  states. 
The  polaron  transformation  applies  a  qubit-state-dependent  shift  to  the  cavity  field,  where  the  complex- valued 
displacement  amplitudes  a j  are  the  cavity  field  amplitudes  conditioned  on  the  qubit  states  and  are  in  general 
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time-dependent.  The  equations  of  motion  for  the  cq  follow  from  the  polaron  transformation  under  the 
constraint  that  the  resulting  Hamiltonian  represent  an  un- driven  cavity.  This  treatment  and  the  resulting 
equations  for  the  cavity  field  are  similar  to  a  semi- classical  treatment  of  the  cavity  dynamics,  while  still  keeping 
track  of  intrinsic  and  qubit-induced  nonlinearities  [52].  The  polaron  frame  best  captures  the  dynamics  of  the  full 
system  when  the  cavity  field  can  be  accurately  described  as  a  coherent  state  whose  complex  amplitude  depends 
on  the  instantaneous  qubit  state.  In  practice  the  use  of  the  polaron  frame  is  restricted  to  weak  measurement 
regime,  as  becomes  evident  when  considering  its  effect  on  the  qubit  ladder  operators: 

=  Dpo_.  (A5) 

Here  the  generalized  displacement  operator  Dp  is  defined  by 

Dp  =  DHa0)D(ai)  =  D^e"^,  (A6) 

with  the  measurement  separability  /?  =  a\  —  a0  and  the  phase  ip  =  Im(ce0af).  This  follows  naturally  from  the 
fact  that  the  polaron  transformation  connects  a  given  qubit  state  with  a  corresponding  coherent  state  in  the 
cavity.  Any  change  in  the  qubit  state  must  be  connected  with  a  simultaneous  change  in  the  cavity  field. 
Importantly,  the  transformation  equation  (A5)  results  in  an  expression  that,  for  non- vanishing  values  of  \/3 1, 
contains  all  orders  of  the  cavity  annihilation  and  creation  operators.  This  makes  it  impractical  to  find  a  reduced 
description  of  the  dynamics  of  the  qubit  alone,  except  for  the  weak  measurement  case,  when  |  f3 1  <C  1  and 
Dp  ~  1.  In  this  case  the  system  reduces  to  a  set  of  equations  for  the  classical  field  amplitudes  cq  conditioned  on 
the  qubit  state  and  a  Hamiltonian  describing  the  qubit  coupled  to  the  olx  and  an  effectively  un-driven  cavity, 
describing  the  quantum  fluctuations  of  the  cavity  field. 

By  applying  the  polaron  transformation  to  the  dissipators,  we  immediately  find  the  measurement- induced 
dephasing  since 


KV[a]p  — >  kD\o\p  +  la 'DWzlp  H — > 

where  Tj  =  \  (3 12  is  the  measurement- induced  dephasing  rate  and  the  omitted  terms  can  be  absorbed  into 

the  coherent  part  of  the  master  equation. 

We  now  find  the  qubit- only  dynamics  by  tracing  the  master  equation  over  the  cavity  degrees  of  freedom  in 
the  Fock  basis  |  n),  to  find  the  time  evolution  of  the  qubit’ s  reduced  density  matrix: 

p  =  J2(n\p\n).  (A7) 

n 

Finally,  the  jump  rate  in  equation  (4)  can  be  found  from  the  qubit’s  Bloch  vector  elements  (crz-).  For  resonant 
qubit  drive,  where  =  £>q,  one  finds 

{<rx)  =  -  (72  +  rd)(<7*), 

(<ry)  =  -  (72  +  rd)(<Ty)  +  Q{a-Z), 

(az)  =  fl(ay)  -  Ti  (<rz)  +  (7f  -  Tl),  (A8) 

where  I]  =  7^  +  7^  includes  all  dissipative  transitions  between  qubit  levels  (i.e.  thermally  induced  by  the 
environment  as  well  as  dressed  dephasing).  Assuming  that  the  measurement- induced  dephasing  dominates, 
such  that  (<jx)  =  (<jy)  =  0,  the  equation  of  motion  for  (  <jz)  takes  the  form: 

(az)  =  -  ( - — —  +  ri)(crz)  +  (7T  -  7y), 

V72  +  rd  ) 

=  -(X\  +  ri)(*z)  +(7t-7|).  (A9) 

We  identify  =  r|5drive  +  T^andTj  =  r^,drive  +  7^  and  find  the  drive -induced  transition  rate  as  in  equation  (4): 

H2 

ff,  drive  f|,  drive  ~~  TTT'  (A  10) 

2(72  +  Id) 


Appendix  B.  Details  of  the  numerical  simulations 

In  contrast  to  the  simple  two-level  qubit  model  used  in  the  previous  description,  the  transmon  qubit  used  in  our 
experiments  is  a  genuine  multi-level  system.  In  the  presence  of  strong  driving  and  measurement,  which  is  the 
regime  of  our  experiments,  the  additional  levels  lead  to  measurably  different  behavior  of  the  coupled  system. 
Our  numerical  simulations  therefore  take  into  account  the  multi-level  nature  of  the  transmon  and  cavity,  using 
experimentally  determined  parameters.  In  the  rotating  wave  approximation,  the  Hamiltonian  is  given  by: 
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H  =  H,+  Ura'a  +  J2sM+a  +  t^V) 

i 

+  ^(6roe~luJ^°tat  +  6*0ela;rofa)  +  ^-(<fqe_la;dtat  +  6*ela;dta),  (Bl) 

where  nw  =  JT  |  i)  (i|  describes  the  transmon  qubit  in  its  eigenbasis  with  eigenstates  |  i)  and  corresponding 

energies  07.  The  qubit  multi-level  ladder  operators  are  defined  as  <7+  =  |i  +  1)  (i|  and  =  cr+ f.  The  cavity 
has  frequency  ujy  and  is  coupled  with  the  coupling  strengths  gj  to  the  qubit  transitions  via  a  Jaynes-Cummings- 
type  interaction.  Each  of  the  g{  can  be  determined  from  the  dipole  moments  corresponding  to  different  qubit 
transitions,  which  are  obtained  by  diagonalizing  the  transmon  qubit  Hamiltonian  in  the  charge  basis  [38] .  The 
cavity  measurement  drive  is  applied  at  frequency  cjro  and  with  amplitude  ero.  In  the  experiment,  driving  of  the 
qubit  is  achieved  as  a  second  order  effect  via  off- resonant  driving  of  the  cavity  at  the  desired  qubit  transition 
frequency;  therefore  in  our  model  the  qubit  is  driven  indirectly  via  the  cavity,  with  drive  frequency  and  drive 

strength  eq.  In  the  dispersive  limit,  (cjr  —  cu10)  =  A  g,  this  term  can  be  written  perturbatively  as  the  directly 
driven  model  OrjZ-  (<7+  +  a® ),  with  the  Rabi  frequencies  for  each  qubit  transition  «  e^gj A. 

However,  we  use  the  exact  model  equation  (Bl)  for  our  simulations. 

To  calculate  the  time  evolution  of  the  system  under  the  influence  of  the  environment,  we  write  the  master 
equation  for  the  density  matrix  as  [30] 

p=  Cp  =  — i [H,  p]  +  nV[a]p  +  2^vT>[U6^]p  +  YjjiV[<T^}]p,  (B2) 

i 

with  photon  loss  from  the  cavity  at  rate  ft,  qubit  decay  from  state  |  i)  at  rate  7-  and  pure  dephasing  acting 
on  the  transmon  due  to  fluctuations  in  its  energy  levels  at  rate  7^.  Here  |i)  (i|  and  = 

x  5  where  T>  is  the  external  magnetic  flux  applied  to  the  qubit  loop.  Then  V  [II^$]  p 

describes  dephasing  of  the  transmon  due  to  slow  fluctuations  in  the  magnetic  field  through  the  qubit  loop.  As 
decay  of  superconducting  qubits  is  mostly  due  to  electromagnetic  coupling  to  spurious  environmental  modes 
[52, 54],  the  relevant  matrix  elements  are  proportional  to  the  respective  dipole  elements  of  the  state  transitions. 
We  therefore  define  the  relaxation  rates  of  higher  transmon  levels  as  7  =  (gj/g0)2Jp  where  ^  is  the  relaxation 
rate  of  the  first  excited  qubit  level.  All  rates  used  in  the  simulations  were  measured  in  independent  experiments. 

To  simulate  the  measurement  process,  we  use  a  stochastic  master  equation  technique  [19, 45, 46].  For  a 
homodyne  measurement,  we  write  the  stochastic  master  equation  as 

P  =  ep  +  Jivrj£(t)(cp  +  ptf  -  (c  +  ct)p),  (B3) 

where  c  =  e^a  is  the  measurement  operator  for  our  system  with  the  local  oscillator  phase  cp  (with  respect  to  the 
cavity  driving  signal),  ft  is  the  measurement  rate  (due  to  leakage  out  of  the  cavity)  and  p  is  the  measurement 
efficiency.  Here  £  (t)  is  a  stochastic  noise  process  defined  by  £  (t)  =  d  W  (t) /dt,  where  d  W  ( t )  is  the  Wiener 
increment  and  dt  is  the  timestep  between  successive  evaluations  of  the  master  equation.  dW  represents  the 
measurement  noise  and  has  the  properties 

(d  W)E  =  0  (d  W2)e  =  dt,  (B4) 

where  the  average  is  over  different  realizations  of  the  noise  process.  We  then  use  the  simulated  stochastic 
dynamics  to  numerically  calculate  the  qubit  populations  ( |  i)  (i|  )t  and  determine  state  transition  rates  with  the 
same  algorithm  used  on  the  experimental  data.  Here  the  average  is  defined  as  (o)t  =Tr{p(f)d}.We  also 
calculate  the  cavity  field  amplitude  (a)t,  the  photon  number  {a^a)t,  and  the  associated  ac  Stark  shifts  to  calibrate 
the  measurement  strength. 

The  numerical  effort  involved  in  these  simulations  is  quite  large,  as  the  strong  measurement  and  qubit  drive 
requires  one  to  take  into  account  up  to  five  levels  of  the  transmon  qubit  as  well  as  up  to  50  photon  states  to  ensure 
proper  convergence.  In  addition,  the  time  steps  used  in  the  stochastic  numerical  simulations  must  be  chosen 
much  smaller  than  any  of  the  intrinsic  timescales  of  the  problem  to  ensure  that  no  unphysical  solutions  will  be 
reached.  Due  to  these  restrictions,  the  number  of  data  traces  from  these  numerical  simulations  is  much  smaller 
than  the  number  of  experimental  traces  to  which  they  are  compared,  and  thus  the  statistical  error  bars  are 
generally  larger. 

Appendix  C.  Multi-level  driven  Purcell  effect 

The  model  Hamiltonian  for  the  Purcell  decay  rate  of  a  multilevel  qubit  in  a  strongly  driven  cavity  is: 

H  =  u)rcfa  +  Vi  |i)(i|  +J2sMia  +  (C1) 

i  i 

where  o[  —  \  i  +  1 )  ( i  |  and  the  sum  goes  over  all  relevant  qubit  levels.  Following  [30, 3 1  ],  we  define  the  qubit’ s 
Purcell  decay  rate  with  the  cavity  in  a  Fock  state  with  photon  number  n  as 
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Figure  7.  Calculated  Purcell  decay  rate  TP  ( n )  for  a  coherent  state  with  mean  photon  number  n,  from  equation  (C3).  The  different 
curves  are  calculated  by  numerical  diagonalization  of  the  Hamiltonian  (Cl)  with  between  two  and  five  qubit  levels  considered.  Plots 
are  made  using  the  same  qubit  and  cavity  parameters  used  in  the  numerical  simulations. 


T{n)  =  n\{g,n\a\e,n)\2,  (C2) 

where  the  eigenstates  \g,  n)  and  |  e,  n)  belong  to  the  subspaces  of  the  Jaynes-Cummings  Hamiltonian  with  n  and 
n  +  1  excitations,  respectively.  If  the  cavity  is  not  in  a  Fock  state  with  definite  photon  number  n,  the  effective 
Purcell  rate  can  be  found  by  averaging  over  the  photon  number  distribution  of  the  cavity  state  P  (n)  as 

r,  =  ^P(n)T(n)/^P(n).  (C3) 


For  a  coherent  state  with  mean  photon  number  n ,  the  photon  number  distribution  is  P  ( n )  =  e~nfin/ n\.  Since  in 
the  absence  of  dissipation  the  Jaynes-Cummings  Hamiltonian  preserves  the  number  of  excitations,  the  task  of 
calculating  the  Purcell  rate  reduces  to  diagonalizing  the  Hamiltonian  for  each  n-excitation  subspace.  The  general 
form  of  the  Hamiltonian  in  the  n-excitation  subspace  can  be  written  as 

'  o  Vng0  0  0 

Vng0  6u>i  Vn  -  Uri  0 

H  -  o  V«  -  lfi  Slo 2  -Jn  -  lg2  ’  ('C4') 

0  0  ■> Jn  —  2g2  8u$ 

\  ...  ...) 


with  the  6c —  kcjr  giving  the  detuning  between  the  kth  qubit  level  and  a  /c-photon-Fock  state.  We  can 
write  the  eigenstates  of  the  n-excitation  subspace  in  the  basis  of  qubit-cavity  product  states  as: 

I g,  n)  =  b^n  | gy  n)  +  b^n  \e,  n  -  1)  +  bin  \ /,  n  -  2)  +•••, 

\e>  n  -  1)  =  cg,n  | g,  n)  +  ce,„  \e,  n  -  1)  +  cf,„  | /,  n  -  2)  +•••,  (C5) 


and  so  on  for  all  other  states  involving  higher  qubit  levels.  Here  |g) ,  |  e),  and  |  /)  denote  the  first  three  qubit 
eigenstates,  and  the  \n  and  c^n  are  complex  numbers  calculated  by  diagonalizing  equation  (C4).  Together  with 
the  action  of  the  photon  annihilation  operator  on  cavity-Fock  states,  a  \  n)  =  sfn  \n  —  1 ),  and  noting  that  a 
will  lower  the  overall  excitation  number  by  one,  we  can  then  identify  the  relevant  matrix  elements  of  a  as 

(g,  n\a  \e>  n)  =  ^Jn  +  1  bgtncgt„+ 1  +  Jn  be,„ce,n+ 1  +  Vn  -  1  bfinCfin+ 1  +•••,  (C6) 

where  we  obtain  the  coefficients  numerically,  using  the  experimental  parameters  for  our  transmon  qubit. 

Figure  7  shows  a  plot  of  the  calculated  Purcell  decay  rates  for  our  experimental  parameters  with  a  coherent  state 
in  the  cavity,  accounting  for  various  numbers  of  qubit  levels.  The  Purcell  decay  is  suppressed  by  the  presence  of 
photons  in  the  cavity,  although  the  effect  is  less  pronounced  when  higher  qubit  levels  are  accounted  for.  The 
calculated  rates  for  more  than  five  qubit  levels  (not  plotted)  are  very  similar  to  those  for  five  qubit  levels. 


Appendix  D.  Maximum  likelihood  expressions  for  rate  estimation 

Our  qubit  is  a  two-state  system  where  state  transitions  obey  Poissonian  statistics,  but  the  finite  bandwidth  of  the 
measurement  apparatus — which  results  from  the  inherent  bandwidth  k  of  circuit  QED  measurement,  the 
bandwidth  of  the  superconducting  parametric  amplifier  and  other  elements  in  the  measurement  chain,  and  the 
Gaussian  filtering  employed  to  increase  SNR  in  post-processing — tends  to  skew  the  observed  dwell  times  in  each 
state  toward  longer  times  and  reduces  the  number  of  events  seen  with  dwell  times  smaller  than  the  measurement 
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A*  B* 


Figure  8.  State  diagram  for  finite  bandwidth  detection.  This  figure  shows  a  system  state  diagram  for  finite  bandwidth  detection.  The 
qubit  is  in  either  state  A  or  B  (circles),  while  the  readout  is  in  state  A*  or  5*.  For  each  qubit/readout  state,  the  likelihood  of 
transitioning  to  a  different  qubit/ readout  state  is  determined  by  the  rates  VAy  VB,  and  Tdet  as  indicated.  Adapted  with  permission 
from  [41]. 


bandwidth.  A  mathematical  model  to  compensate  for  these  effects  has  been  presented  by  Naaman  and 
Aumentado  [43]. 

Their  model  allows  the  states  of  the  qubit  (denoted  A  and  B )  to  be  independent  of  the  corresponding  states  of 
the  readout  (denoted  A*  and  £*).  A  diagram  of  this  model,  adapted  from  [43],  is  shown  in  figure  8.  For  each 
readout  state,  the  qubit  can  be  in  either  state  A  or  B.  Transitions  occur  between  A  and  B  with  rates  VA  and  VB. 
However,  when  the  qubit  and  readout  ‘disagree’,  i.e.  states  ( B ,  A*)  and  (A,  £*),  the  readout  can  also  make  a 
transition  (to  (B,  5*)  or  (A,  A*),  respectively)  with  rate  VBet .  Assuming  that  the  readout  does  not  change  states 
unless  the  qubit  has  changed  states  (no  false  positives),  one  can  write  the  probability  distribution  function  h  ( t ) 
for  the  observed  dwell  times  in  A*  as  a  function  of  the  underlying  rates  Ifo  Tg,  and  [43]: 

Hr,  rA,  rB,  rdet)  =  TrArdete-At/2sinh^),  (di) 


where  A  =  VA  +  Tg  +  F^et  and  0A  =  ^  A2  —  4TATBet .  The  probability  distribution  for  dwell  times  in  state  B*  is 
identical  but  with  VA  and  VB  interchanged  (defining  0B  =  >/A2  —  41^1^  )>  as  can  be  seen  from  the  symmetry  of 
the  model.  We  can  calculate  the  complementary  cumulative  distribution  function,  often  called  the  survival 
function,  by  integrating  the  expression  in  (Dl): 


s(t;  Ia,  Tg,  Tjet) 


i: 


h(t')dt'  = 


p-Af/2 


h  sinh  +  H  cosh  ^  j 


(D2) 


The  survival  function  is  useful  for  handling  right- censored  dwell  times,  where  the  observed  dwell  time 
represents  a  lower  bound  on  the  true  dwell  time  in  the  state.  The  survival  function  s  (r )  gives  the  probability  of  a 
dwell  time  t  >  r,  which  is  precisely  the  probability  distribution  function  needed  for  right- censored  events. 
Given  a  data  set  of  nA*  dwell  times  in  state  A*  and  nB *  dwell  times  in  state  B*,  denoted  { t{\  and  { ty},  respectively, 
as  well  as  corresponding  censoring  variables  { <%}  and  { <5;  },  where  6  =  0(1)  indicates  an  uncensored  (right- 
censored)  dwell  time,  we  can  write  the  likelihood  function  as: 

nA* 

£(Ta,  Tg,  Tdet)  =  Yi  hitil  Ia,  Tg,  Tdet)1-^  s(til  Ta,  Tg,  Tdet)^' 

i=  1 

nB* 

x  n  h(tp  rB»  rd,  rB,  rA,  rdet)^.  (D3) 

;=  i 


We  then  use  nonlinear  optimization  methods  to  maximize  the  value  of  C  by  varying  the  parameters  VAy  Tg, 
and  Tjet  •  Because  of  the  limitations  of  floating  point  arithmetic,  we  actually  perform  the  maximization  on  the 
log-likelihood  function  L  =  In  (£),  which  has  the  form: 

nA* 

L( Ta,  Tg ,  Tdet)  =  X][(l  —  $i) In [h (tz-;  VAi  T5,  Tdet)]  +  <5zTn[ste  TA,  VB,  Idet)] ] 

i=  1 

nB * 

+  X3[(1  -  A  in  [h(tp  rB,  rA,  rdet)]  +  s,  in  [s(tr,  rB,  rA,  rdet)]].  (D4) 

7=1 

The  functional  forms  of  In  [h  (t)  ]  and  In  [s  ( t )  ]  for  dwell  times  in  state  A*  are: 

ln[h(P,  rA,  rB,  rdet)]  =  In |  j  -  y  +  ln[sinh(y“)j’  (D5) 
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In  [s(t;  Ia,  Tb,  Tjet)] 


—In  (0A)  —  +  lnj^A  sinh^-^j  +  0A  cosh 


(D6) 


The  corresponding  functions  for  dwell  times  in  state  5*  can  be  found  by  interchanging  VA  with  Tg,  and  0A 
with  0B.  For  sufficiently  large  values  of  0A  t,  these  expressions  may  cause  overflows  in  double-precision  floating 
point  arithmetic.  Therefore,  for  0At  >  40,  we  make  the  following  approximations,  which  introduce  fractional 
errors  of  less  than  10-18  for  each  approximated  term: 


lnj^Asinh^-^-j  +  #ACOsh|-^-j 


+  In  (A  +  9a)  —  In  2. 


(D7) 

(D8) 
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